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ABSTRACT 


Waves  in  the  surf  zone  are  a highly  nonlinear  process 
which  is  evident  by  the  appearance  of  secondary  waves.  The 
secondary  waves  appear  as  strong  peaks  in  the  period  PDFs 
corresponding  to  the  first  harmonic  of  the  peak  of  the  wave 
spectrum.  The  strong  first  harmonic  period  peak  is  also 
reflected  in  the  highly  correlated  height  and  velocity  PDFs. 

Ereaking  waves  may  be  viewed  as  a gradation  between  two 
extremes.  The  period,  height  and  velocity  PDFs  for  spilling 
breakers  generally  were  unimodal,  whereas  for  plunging 
breakers  strong  bimodality  was  found. 

The  joint  probability  density  functions  for  periods  and 
heights  of  the  breaking  waves  show  high  correlation  (0.60- 
0.30)  which  says  that  greater  wave  periods  are  associated 
with  larger  breaker  heights.  The  joint  PDFs  of  period  and 
particle  velocity,  and  velocity  and  height,  suggest  that  the 
maximum  onshore  particle  velocities  are  correlated  with  both 
the  wave  periods  and  wave  heights. 

Due  to  the  high  probability  of  the  secondary  waves,  the 
mean  wave  period  for  breakers  is  a poor  descriptor  of  the 
average  period  of  the  offshore  incident  waves. 
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Waves  in  deep  water  can  generally  be  represented  to  a 
first  approximation  as  a Gaussian  process  having  zero  skew- 
ness and  kurtosis  equal  to  three.  The  random  sea  surface 
is  viewed  as  the  linear  superposition  of  statistically  inde- 
pendent spectral  components  which  is  adequately  described 
by  the  energy  density  spectrum.  The  restriction  of  statis- 
tical independence  implies  the  wave  components  act  as  free 
waves  and  there  is  no  coupling  or  interaction  between 
various  spectral  components. 

As  the  waves  begin  to  shoal,  the  crests  become  more 
peaked,  and  the  troughs  become  more  elongated;  this  results 
in  a positive  value  of  the  third  statistical  moment,  i.e., 
the  probability  density  distribution  of  sea  surface  eleva- 
tions becomes  skewed  to  the  right.  The  skewness  is  a result 
of  phase  couplirfg  of  the  wave  components,  which  is  indica- 
tive of  the  increasing  importance  of  the  non-linearities 
of  the  system.  This  initially  weakly  non-linear  system 
found  in  shallow  water  may  be  adequately  described  by  a 
Gram-Charlier  probability  distribution  (Longuet-Higgins , 
1963) . 

Near  breaking  waves  become  asymmetrical  about  the  crest 
with  a continually  steepening  of  the  wave  face.  Miller 
(1976)  has  shown  that  a jet  is  then  formed  at  the  wave 
crest  at  breaking,  which  extends  forward  and  then  down 
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until  it  closes  with  the  forward  face  of  the  breaking  wave 
form.  The  asymmetry  about  the  wave  crest  yields  a change 
in  the  fourth  moment,  or  kurtosis,  which  is  also  indicative 
of  the  increasing  non-linearities.  Thornton  et  al.  (1976) 
have  noted  that  secondary  wave  formation  becomes  apparent 
on  the  trailing  face  of  the  wave  form,  due  to  transfer  of 
energy  from  the  primary  frequency  to  its  harmonics  which 
also  points  to  a strongly  nonlinear  system. 

This  thesis  uses  data  collected  as  described  by  Galvin 
(1975)  and  Olsen  (1977)  plus  additional  data  collected  by 
the  author.  The  experiments  were  designed  to  measure  the 
sea  surface  elevation  and  the  onshore  particle  velocities 
by  means  of  capacitance  wave  gauges  and  electromagnetic 
flowmeters,  respectively.  The  experiments  were  conducted 
at  a variety  of  beaches  in  order  to  include  wave  types  of 
spilling,  plunging,  and  collapsing  forms.  The  characteris- 
tics of  the  breaking  wave  forms  are  a function  of  the  off- 
shore wave  spectrum  and  the  nearshore  bottom  slope  and 
topography,  where  the  increasing  beach  slope  increases  the 
rapidity  of  the  breaking  of  the  waves. 

Empirical  distributions  were  obtained  for  the  periods, 
wave  heights,  and  maximum  onshore  particle  velocities,  for 
waves  which  had  been  defined  by  the  zero-up-cross  method. 
Each  of  these  distributions  was  compared  with  the  Rayleigh 
probability  density  function  (PDF).  Joint  PDFs  of  periods 
and  wave  heights,  maximum  velocities  and  wave  heights,  and 
periods  and  maximum  velocities  were  calculated. 
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II.  HISTORICAL  BACKGROUND 

A.  QUALITATIVE  DESCRIPTION 

Thornton  et  al.  (1976)  have  described  a number  of  ob- 
served universal  similarities  for  the  various  types  of 
breakers  occurring  on  different  beaches.  Normally  there  is 
a quick  drawdown  of  water  just  before  the  breaker  arrives, 
followed  by  a steep,  nearly  vertical  leading  edge.  Miller 
(1976)  suggests  that  subsequent  to  the  development  of  the 
leading  vertical  edge,  a jet  forms  at  the  crest.  This  jet 
progressively  extends  itself  forward  and  down  until  it  unites 
with  the  forward  slope  of  the  breaking  wave,  resulting  in 
the  formation  of  a vortex  with  a circulation  around  a 
cavity.  When  the  tip  of  the  "jet"  touches  the  forward  slope 
of  the  breaking  wave,  the  resulting  force  forms  a second 
weaker  vortex  in  front  of  the  first  one.  In  this  manner,  a 
series  of  vortices  yields  a series  of  large,  single  pulse 
eddies,  which  contribute  to  the  rapid  decay  of  the  wave 
during  breaking.  It  has  been  found  that  the  magnitudes  of 
breaker  vortices  are  generally  correlated  with  a gradation 
from  plunging  to  spilling  breakers. 

Thornton  et  al.  (1976)  have  observed  a sloping  profile 
toward  the  trailing  edge  of  the  wave.  Secondary  waves  are 
frequently  noted  on  the  trailing  edge  both  on  the  analog 
strip  charts,  Figure  1,  and  by  visual  observations  in  the 
field.  The  secondary  waves  have  been  observed  to  grow  during 
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the  shoaling  process,  especially  just  prior  to  breaking 
when  the  primary  waves  are  steepening  significantly.  The 


secondary  waves  are  a result  of  the  transfer  of  energy  from 
the  primary  wave  frequency  to  harmonics,  and  help  maintain 
the  potential  energy  across  the  surf  zone. 

Gallagher  (1971)  has  photographed  secondary  wave  forma- 
tion occurring  across  reefs.  He  suggests  that  the  secondary- 
waves  are  brought  about  by  a response  to  vertical  accelera- 
tion near  the  front  of  a broken  or  breaking  swell.  When  the 
amplitude  and  surface  curvature  of  a shallow  water  wave  ex- 
ceed certain  values,  the  non-linearities  of  the  vertical 
accelerations  bring  soliton  wave  forms  into  being. 

Figure  1 shows  that  the  maximum  particle  velocity  occurs 
near  the  time  of  maximum  surface  elevation.  It  can  also  be 
noted  in  Figure  1 that  there  are  secondary  velocity  maxima 
which  correlate  with  the  secondary  surface  elevation  maxima. 
The  secondary  velocity  peaks  also  are  a result  of  the  non- 
linearity associated  with  the  breaking  wave. 

B.  ENERGY  DENSITY  SPECTRA  OF  BREAKING  WAVES 

In  order  to  better  understand  the  PDFs  as  determined  in 

this  study,  the  energy  density  spectra  measured  by  Galvin 

(1975)  and  Olsen  (1977)  were  reviewed.  The  energy  density 

spectrum  is  proportional  to  potential  wave  energy  within  a 

frequency  band.  The  energy  density  spectrum,  G (f) , indicates 

2 

how  the  variance,  ax  , is  distributed  over  frequencies,  f, 
where 
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CO 


cx(f)df 


(1) 


(Kinsman,  1965), 

The  energy  density  spectra,  as  estimated  from  the  sea 
surface  elevation  records,  were  used  to  find  the  frequency 
of  peak  density  and  the  occurrence  of  harmonics.  The  power 
spectra  that  were  measured  by  Galvin  showed  a peak  density 
over  the  frequency  range  of  0.061  to  0.092  Hz  (16.4  to  10.9 
sec).  First,  second  and  third  harmonics  can  be  observed 
as  shown  in  the  spectra  of  surface  elevation  and  horizontal 
velocities  for  4 March  1975.  Additionally,  subharmonics 
were  observed  at  0.011  and  0.922  Hz  in  the  velocity  spectra. 
Galvin  (1975)  has  suggested  that  the  observed  low  frequency 
sub-harmonic  peaks  may  be  due  to  edge  waves  in  the  surf  zone. 

Thornton  et  al.  (1976)  note  that  the  appearance  of  strong 
harmonics  could  be  real  due  to  secondary  waves  or  artificial 
due  to  Fourier  analysis  of  a peaked  wave  form.  It  is  noted 
that  secondary  maxima  of  n are  apparent  in  the  strip  chart 
record  shown  in  Figure  1.  It  is  felt  that  the  secondary 
spectral  peaks  are  real  since  secondary  waves  appear  on  the 
strip  charts. 

C.  PROBABILITY  DENSITY  FUNCTIONS  OF  PERIODS  AND  HEIGHTS 

Studies  of  the  PDFs  of  period,  height,  and  wave  lengths 
have  been  conducted  over  the  last  20  years.  The  majority 


u 
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Figure  2.  Energy  Density  Spectra  of  n and  v,  4 March  1975. 
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of  these  studies  have  been  mainly  concerned  with  theoretical 
or  deep  water  aspects  of  the  problem. 

The  theoretical  distribution  of  wave  amplitudes  for  a 
narrow-banded  Gaussian  surface  elevation  in  deep  water,  was 
shown  by  Longuet-Higgins  (1952)  to  be  a Rayleigh  distribu- 
tion. Tayfun  (1977),  in  studying  the  transformation  of  deep 
water  waves  to  shallow  water  waves,  showed  that  the  Rayleigh 
distribution  for  wave  amplitude  was  generally  applicable  to 
all  bandwidths.  Using  the  assumption  that  the  wave  height 
is  twice  the  wave  amplitude,  the  wave  height  PDF  is  then  re- 
presented by: 

r t 

p (II) dH  = ^ exp  j^l  dll  (2) 

where  a is  the  variance  of  the  surface  elevation.  Using 
pressure  records  in  the  Gulf  of  Mexico,  Longuet-Higgins 
(1975)  observed  that  the  R.ayleigh  distribution  fit  the  ob- 
served distribution  reasonably  well.  He  found  that  there 
is  a slight  excess  of  waves  with  heights  near  the  middle 
of  the  range  and  a deficit  at  the  two  extremes.  Since  much 
of  the  high  frequency  portions  of  the  wave  records  are  fil- 
tered out  by  the  pressure  transducer,  Longuet-Higgins  (1975) 
suggests  that  the  narrow-band  approximation  may  not  be  as 
applicable  for  the  unfiltered  records.  In  shallow  water 
with  much  steeper  waves,  the  distribution  can  again  be 
expected  to  be  less  applicable  due  to  the  non-linearities 
which  become  increasingly  important. 

15 


A 


Since  Che  statistics  of  H and 


for  a broad-banded 


spectrum  were  still  somewhat  unclear,  Goda  (1970)  simulated 
wave  profiles  using  an  infinite  series  where  the  amplitudes 
were  specified  by  various  theoretical  spectra  and  the  phase 
was  random  with  a computer.  He  then  examined  the  simulated 
wave  records  for  surface  e levation , cres t- to- trough  wave 
heights  and  zero-up-cross  wave  heights.  He  found  that  for 
the  zero-up-crossing  determination  of  wave  heights,  the  Ray- 
leigh distribution  is  a good  approximation  irrespective  of 
the  spectral  width. 

Goda  (1970)  suggests  that  the  deviation  of  the  wave 
height  distribution  from  the  Rayleigh  type  in  shallow  water 
should  be  considered  separately  from  the  change  in  the  func- 
tional shape  of  the  wave  spectrum.  It  must  be  a highly  non- 
linear phenomenon,  which  requires  full  analysis  with  the 
theory  of  non-linear  wave  spectra. 

Bretschneider  (1959)  investigated  the  PDFs  of  wave 
heights,  H,  wave  lengths,  L,  and  wave  periods,  T.  He  con- 
cluded from  his  empirical  study  that  in  deep  water  wave 
systems,  the  PDFs  of  wave  heights  and  lengths  can  be  approxi- 
mated by  a Rayleigh  distribution.  Using  the  relationship, 
wave  length  proportional  to  period  squared  for  deep  water 
linear  waves,  he  concluded  that  the  PDF  of  periods  squared 
is  a Pv-ayleigh  PDF  given  by 


x2 


r— 

rTl~| 

exp 

-0.675 

tJ 

(3) 


where  T is  the  mean  period. 
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Using  simulated  wave  data,  Goda  (1970)  found  that  the 
ores t- to- ores t periods  squared  were  generally  not  well  ap- 
proximated by  the  Rayleigh  PDF.  The  mean  and  standard  devia- 
tions of  the  periods  by  this  method  were  found  to  be  a 
function  of  the  spectral  width.  He  found  that  the  periods 
determined  by  the  zero-up-cross  method  are  only  in  fair 
agreement  with  Bretschneider  (1959).  If  the  waves  are  of 
the  type  which  have  sharply  peaked  spectra,  double  peaked 
spectra,  or  a flat  spectrum,  Bretschneider' s (1959)  distri- 
bution was  totally  in  disagreement  with  that  of  the  simu- 
lated waves . 

Contrary  to  Bretschneider  (1959)  and  Goda  (1970), 
Longuet-Higgins  (1975)  noted  that  in  general  the  deep  water 
wave  periods  were  bell  shaped  (but  not  Gaussian).  He  defines 
the  period  PDF  when  calculated  by  the  zero-up-cross  method 
as  , 

p(n)  = j ; n = (t-<t>)/v<t>  (4) 

2 (1+n  ) Z 

where  <t>  is  the  mean  period  and  v is  proportional  to 
the  spectral  width  where 

v = (u2/u0)^C<t>/2^)  (5) 

and  uQ  and  u2  are  c^e  zero  and  second  moments  of  the 
energy  spectrum. 

Koele  and  de  Bruyn  (1964),  Goda  (1967),  and  Siefert 
(1970)  have  observed  that  the  distribution  of  wa^e  heights 
in  shallow  water  does  not  correlate  well  with  the  Rayleigh 
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distribution.  Kuo  and  Kuo  (1974)  have  suggested  that  this 
is  due  to:  1)  the  non-linear  effects  of  wave  interactions 
yielding  more  larger  waves,  2)  the  effect  of  bottom  friction 
yielding  a reduction  in  the  low  frequency  components,  and 
3)  the  effects  of  wave  breaking,  which  would  truncate  the 
distribution  and  transfer  some  of  the  kinetic  energy  to  the 
high  frequency  components.  Therefore  the  PDFs  of  wave 
heights  with  a certain  intensity,  may  be  considered  to  ap- 
proximate a conditional  Rayleigh  distribution  truncated  by 
the  local  breaking  height. 

A zero-up- cross  analysis  was  conducted  on  waves  measured 
on  a reef  by  Black  (1978),  who  then  compared  the  observed 
distribution  with  the  Rayleigh,  truncated  Rayleigh,  and 
the  Ueibull  distributions.  He  observed  that  the  truncated 
Rayleigh  distribution  fit  well  for  those  waves  in  the  breaker 
zone  when  the  heights  were  depth  limited.  The  truncated 
Rayleigh  distribution  did  not  fit  well  for  those  waves 
measured  either  inshore  or  offshore  of  the  breaker  zone. 

He  has  found  that  a Ueibull  distribution  more  closely  approx- 
imated these  measured  distributions. 

It  would  be  expected  that  the  Ueibull  distribution  would 
yield  the  better  fit.  The  Rayleigh  distribution  is  a func- 
tion of  the  variance  of  the  data,  whereas  the  Ueibull  dis- 
tribution is  a function  of  higher  moments  about  the  mean. 
Neither  theoretical  distribution  adequately  describes  ob- 
served distributions  which  are  multi-modal. 

Black's  (1973)  observations  of  wave  periods  showed  a 
distribution  which  was  skewed  to  the  right.  This  does  not 
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fit  well  with  the  symmetrical  Longuet-Higgins  (1975)  distri- 
bution. It  was  found  that  the  period  distribution  more 
closely  followed  the  Rayleigh  distribution  than  the  period 
square  as  proposed  by  Bretschneider  (1959).  It  was  found 
that  the  Ueibull  distribution  most  closely  fit  the  empiri- 
cal PDF,  but  there  is  considerable  variation  in  the  peaked- 
ness parameter.  His  attempts  in  relating  the  period  statis- 
tics to  the  position  on  the  reef  were  unsuccessful.  He 
feels  that  this  is  due  to  the  non-linearities  of  the  waves 
in  shallow  water. 

In  this  thesis,  the  empirical  data  collected  for  breaking 
waves  were  compared  with  the  Rayleigh  distribution. 

D.  JOINT  PROBABILITY  DENSITY  FUNCTIONS 

A study  of  the  joint  PDF  for  wave  heights,  ",  and  lengths, 
X,  which  can  be  viewed  as  heights  and  periods)  was  accom- 
plished by  Bretschneider  (1959).  The  general  form  of  the 
joint  PDF  for  two  related  variables  (n,X)  is  the  product  of 

the  marginal  PDF  (p(n))  and  the  conditional  PDF  (p(\)). 

n 

p (n,  X)  = p (n)  •pr|(X)  (6) 


then 


P (n) 


00 


(n, X) dx 


o 


( X ) dX 


(7) 


Assuming  that  the  marginal  PDFs  are  Rayleigh,  Bret- 
schneider (1959)  examined  the  extreme  cases  of  0 and  +1  cor- 
relation. For  the  case  of  zero  correlation,  the  marginal 
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Rayleigh  PDFs  are  given  by 


P ( 

p)  = t n 

-7 

exp 

: n = H/H 

(8) 

P ( 

rr 

')  = 2 X 

2 

exp(-iTX  / 4) 

; X = L/D 

(9) 

For 

zero 

corre lation 

.,  i.e.  total 

independence  (Figure 

3)  , 

the 

j oin 

t PDF  is. 

P (n, x) 

= p (n)  *p  (X) 

(10) 

so 

that 

the 

joint  PDF  of  heights 

and  lengths  is  given 

as, 

p (n, 

X)  ■ 

rr2 

= 4 n exp 

(-irn2/4)  • X e 

xp  (-tX2/4) 

(ID 

In 

terms 

of 

heights 

and  periods  , 

equation  15  is , 

P (n , 

T)  - 

= 1.35  Te 

2 

Xp  (-  TT|  /4)  • 

t3  exp (-0 . 6 75  t4) 

(12) 

where  t = T/T. 

For  the  case  of  total  dependence,  the  correlation  coef- 
ficient is  equal  to  one.  For  a correlation  coefficient  of 
one,  all  data  points  on  a plot  of  joint  Rayleigh  PDFs  fall 
on  a 45  degree  straight  line  passing  through  the  origin, 
n=-\=0 . Assuming  that  a joint  Rayleigh  distribution  is  appl 
cable,  the  following  relation  applies: 

p(n,X)  = j n exp(-:rn2/4)  • £ X exp(--TX2/4)  (13) 

which  is  plotted  in  Figure  3. 
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Figure  3.  Statistically  Independent  Joint  Rayleigh  PDF. 


Collins  (1970)  studied  the  PDFs  of  breaking  waves  de- 
rived from  deep  water  waves  having  a joint  Rayleigh  distri- 
bution of  heights,  lengths  and  angles.  He  assumed  that 
H = yd  , where  II  is  the  breaker  height,  d is  the  breaker 
depth,  and  y is  the  breaking  criterion  assumed.  The  limit- 
ing factor  of  waves  breaking  in  the  surf  zone  yields  a dis- 
tribution better  approximated  by  a truncated  Rayleigh  dis- 
tribution. As  expected,  he  found  that  the  longest  deep  water 
waves  were  associated  with  the  largest  breaking  waves. 

Collins  (1970)  concluded  from  the  study  of  the  joint 
PDFs  of  wave  heights  and  lengths,  that  the  effect  of  deep 
water  angle  of  approach  is  relatively  minor  for  breaking 
waves.  The  breaking  wave  steepness  was  found  to  be  very 
sensitive  to  the  choice  of  breaking  criteria. 

A theoretical  expression  for  the  joint  distribution  of 
wave  periods  and  amplitudes  was  derived  by  Longuet-Higgins 
(1975)  for  a narrow  frequency  spectrum 

2 

p(C,n)  = — 2 — n exp  (-C“  (1+0/ 2)  (14) 

(2r)'2 


where  5 and  n are  the  normalized  amplitudes  and  periods 
as  defined  by, 


5 = 


(15) 


and 


n = (t-<t>)/v<t> 


(15) 


where  \>  is  defined  in  equation  5.  The  wave  amplitude 


is  a , y is  the  zero  moment  of  the  energy  spectrum, 
and  <t>  is  the  mean  wave  period.  Longuet-Higgins  (1975) 
compared  this  theoretical  distribution  of  ocean  waves  to 
that  of  the  deep  water  data  obtained  by  3re ts chneider  (1959) 
and  found  good  agreement. 

Chakrabarti  and  Cooley  (1977)  compared  the  theoretical 
joint  distribution  of  Longuet-Higgins  (1975)  with  wave  data 
recorded  during  a 1961  storm  in  the  North  Atlantic.  The 
spectrum  of  the.e  waves  did  not  fall  in  the  category  of 
narrow-band.  Through  the  comparison,  it  was  found  that 
there  was  a considerable  agreement  at  the  higher  wave  heights, 
but  a definite  trend  away  from  the  theoretical  curve  at  the 
lower  wave  heights. 
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III.  EXPERIMENT 


A.  EXPERIMENTAL  SITES 

Three  different  California  experimental  sites  were  used 
to  obtain  data  on  the  three  types  of  breaking  waves--plung- 
ing,  spilling,  and  collapsing.  During  the  period  4 March  to 
10  March  1975,  plunging  and  spilling  waves  were  measured  by 
Galvin  (1975)  on  Del  Monte  3each,  Monterey,  California. 

This  locale  was  chosen  because  the  waves  are  highly  refracted 
and  directionally  filtered,  yielding  narrowbanded  swell  off- 
shore which  impinges  almost  perpendicularly  on  the  beach. 

On  29  May  1975  Galvin  (1975)  conducted  a second  experi- 
ment on  the  Carmel  River  Beach,  which  is  located  approximate- 
ly five  miles  south  of  Monterey,  California.  This  area  is 
in  an  embaynent , where  the  waves  are  again  highly  refracted 
and  directionally  filtered.  The  beach  at  this  site  is  very 
steep  resulting  in  plunging  and  often  collapsing  and  surging 
breakers.  Due  to  the  steepness  of  the  beach  and  the  rapidity 
of  breaking  of  the  shoaling  waves,  reflected  waves  are  found 
to  be  present. 

Spilling  breakers  on  a gently  sloping  beach  were  measured 
in  March  1977  at  Torrey  Pines  Beach,  a site  just  north  of 
La  Jolla,  California  (Olsen,  1977). 

On  3 and  9 March  1978  data  were  collected  by  this  author 
on  Del  Monte  Beach,  Monterey,  California.  On  the  first  day 
of  the  experiment,  very  narrowbanded  swell  was  observed 
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breaking  on  the  shore,  mainly  of  the  plunging  type.  A storm 
arrived  late  on  3 March  and  continued  through  9 March.  Winds 
of  approximately  15  knots  yielded  a very  confused  sea. 
Breakers  were  of  both  the  spilling  and  plunging  type  and 
were  observed  to  break  at  various  small  angles  relative  to 
the  beach.  Sea  surface  elevations  and  velocities  were  col- 
lected for  approximately  three  hours  on  each  of  these  days 
on  both  the  flowing  and  ebbing  sides  of  high  tide. 

B.  INSTRUMENTS 

Capacitance  type  gauges  were  used  to  measure  the  change 
in  the  sea  surface  elevation  as  explained  by  Thornton  et  al. 
(1976).  The  wave  gauges  were  constructed  of  3/3  inch  poly- 
propylene covered  stainless  steel  rod.  Calibration  was 
accomplished  in  the  laboratory/  statically  by  lowering  the 
wave  gauge  into  a deep  water  tank.  Depth  of  immersion  was 
then  plotted  against  volts. 

The  velocity  of  the  water  particles  was  measured  with 
Marsh-McBimey  Model  721  and  722  Electromagnetic  Current 
meters.  These  flowmeters  work  on  Faraday's  principle  of 
electromagnetic  induction  and  are  capable  of  recording  the 
particle  velocities  in  orthogonal  directions  in  a plane. 

The  meters  were  dynamically  calibrated  in  the  laboratory. 
Calibration  was  accomplished  by  oscillating  the  current 
meters  in  a water  tank  with  a prescribed  motion  simulating 
the  flow  conditions  which  would  be  expected  to  be  found  in 
the  surf  zone.  The  flowmerer  voltage  output  was  found  to 
be  linear  in  the  velocity  ranges  found  in  the  breakers. 


25 


All  instruments  were  placed  on  towers  in  the  surf  zone 
at  low  tide  at  the  expected  point  of  breaking  of  the  waves 
for  the  subsequent  high  tide.  Each  tower  held  a capacitance 
wave  gauge  and  an  electromagnetic  flowmeter.  The  flowmeter 
was  placed  below  the  wave  gauge,  approximately  0.5  m from 
the  bottom,  to  insure  constant  submergence.  Data  were  col- 
lected within  an  hour  or  two  on  both  the  ebb  and  flow  sides 
of  high  tide.  Analog  recordings  were  made  on  both  a strip 
chart  and  magnetic  tape. 

C.  DATA  ANALYSIS 

All  signals  were  digitized  from  their  analog  form  by 
means  of  a Vidar  Corporation  32-channel  digital  data  acqui- 
sition system.  The  digitized  data  were  then  processed  on 
an  1 31 1 360/67  computer.  Reco.rd  lengths  of  approximately 
30  minutes  from  each  data  set  were  analyzed.  Specific  de- 
tails of  the  analysis  can  be  found  in  Appendix  A. 

The  data  were  edited  by  excluding  any  value  greater 
than  six  standard  deviations  from  the  mean  and  replacing 
that  value  with  the  previous  value.  The  values  were  then 
linearly  detrended  to  exclude  the  effects  of  the  rising  or 
falling  tides.  The  mean  and  variance  of  sea  surface  eleva- 
tions were  calculated  for  each  record  as  well  as  the  mean 
and  variance  of  the  heights , periods  and  maximum  velocities 
all  of  which  are  defined  on  page  27. 

It  was  desired  to  only  examine  the  sea-swell  band  of 
frequencies  lying  between  0.03-1.0  Hz  (30  sec-1  sec).  Hence, 
the  data  were  passed  through  both  a high  and  low  pass  digital 
filter. 
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The  high  pass  filter  had  a cut-off  frequency  of  0.03  Hz 
(30  sec).  The  high  pass  filter  used  a Fast  Fourier  Trans- 
form (FFT)  algorithm  to  obtain  the  Fourier  spectrum  of  the 
entire  30  minute  record.  The  Fourier  coefficients  correspond- 
ing to  0-0.03  Hz  were  then  used  to  synthesize  a low  frequency 
time  series  which  was  subtracted  from  the  wave  record. 

The  low  pass  filter  was  a 25  weight  inverse  transform 
filter  (Davidson,  1970),  as  shown  in  Figure  4.  A cut-off 
frequency  of  1.2  Hz  was  used  with  a terminal  frequency  of 
1.3  Hz. 

Wave  heights  were  determined  from  the  surface  elevation 
record  by  means  of  the  zero-up-cross  method.  The  zero-up- 
cross  technique  defines  a wave  period  as  the  interval  between 
adjacent  upcrosses.  The  maximum  and  minimum  of  the  surface 
elevation  within  the  interval  defines  the  crest  and  trough 
of  a wave.  Maximum  onshore  particle  velocities  were  calcu- 
lated from  the  flowmeter  records.  Maximum  onshore  velocity 
was  defined  as  the  maximum  velocity  amplitude  occurring  be- 
tween successive  zero  upcrossings  as  defined  by  the  sea 
surface  elevation  record.  Joint  PDFs  of  wave  height  and 
period,  wave  height  and  maximum  velocity,  and  maximum  veloc- 
ity and  period  were  calculated  from  the  analyses. 

Heights  and  maximum  velocities  were  normalized  by  divid- 
ing by  their  respective  standard  deviations.  The  periods 
were  normalized  by  dividing  by  their  mean  period.  The  nota- 
tion us e d is 


H = H' /a 

n 

V = V'/o 

v 

T = T' /? 


(17) 


FREQUENCY  (HZ. 


Figure  4. 


Low  Pass  Digital  Filter. 
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where  heights  (H1 ) , maximum  velocities  (V' ) , and  periods 
(T ' ) are  the  values  as  determined  by  the  zero-up-cross 

method.  f is  the  mean  period,  and  a and  a are  the 

n v 

standard  deviations  of  the  surface  elevations  and  particle 
velocities  respectively.  Histograms  were  calculated  for 
the  normalized  heights,  periods,  and  maximum  velocities. 
The  histograms  were  transformed  into  PDFs  by  dividing  the 
number  of  value  over  an  interval  by  the  total  number  of 
values  in  the  array  and  by  dividing  by  the  interval  width. 
The  mathematical  expression  for  the  exact  PDF  is 

p(x)  - lim  0 1 (Um  .Px|,  (18) 

Ax  [N  j v ' 


where  is  the  number  of  times  the  variable  occupies  an 

interval,  N is  the  total  number  in  the  record  observed, 

Ax  is  the  interval  width  and  p(x)  is  the  PDF.  A restric- 
tion on  the  probability  density  is  that 

-CO 

J ? (x)  dx  = 1 . 

o 

Joint  probability  density  functions  of  periods  versus 
heights,  velocities  versus  heights  and  periods  versus 
velocities  were  calculated  in  a similar  fashion.  A 15  by  15 
grid  was  used  to  determine  the  frequency  of  values  which 
simultaneously  fill  each  joint  interval.  Mathematically 
the  probability  that  the  arrays  x and  y simultaneously 
assume  values  within  the  range  (x,  x+Ax) , (y,  v+Ay)  may  be 


obtained  by  taking  the  ratio  of  the  total  number  of  times  of 
each  occurrence  within  a bin,  N , and  the  total  number  of 
observations,  N . The  exact  probability  density  is  approached 
as  the  number  of  observations  approaches  infinity  the  interval 
widths  approach  zero.  The  joint  PDF  can  be  defined  as: 


P (x,y) 


lira  1 rlim  I'  xy  [ - 1 

x 00  AxAy  N -►  30  [N  J ' 

y ~ 30 


(19) 


with  the  restriction 

OO  +QO 

J p(x,y)  dx  dy  = 1 

o o 


The  energy  density  spectrum  (power  spectrum)  was  obtained 
by  means  of  a Fast  Fourier  Transform  (FFT)  algorithm.  The 
energy  density  spectrum  tells  how  the  variance  is  distributed 
with  respect  to  frequency. 

The  spectral  width  parameter  was  calculated  using  the 
equation : 


lA0rilU~m2 

rn0ni4 


(20) 


where  , m2 , and  are  the  spectral  moments  (Cartwright 

and  Longuet-Higgins , 1956).  For  s=0  the  waves  have  only  one 
frequency  and  the  energy  spectrum  would  be  represented  by  a 
spike.  As  £-1  the  sea  surface  can  be  described  as  a broad- 
band process  for  which  the  distribution  would  be  Gaussian 
and  the  spectrum  would  approach  a horizontal  line.  Swell 
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can  therefore  be  described  as  a narrow-band  process  and  seas 
as  broad-band.  Due  to  the  appearance  of  secondary  waves,  a 
narrow-banded  deepwater  spectrum  of  waves  broadens  in  the 
surf  region. 

The  correlation  coefficients  of  the  joint  values  were 
calculated  to  obtain  a measure  of  the  amount  of  linear  rela- 
tionship of  random  pairs  of  values  of  periods  and  heights, 
velocities  and  heights,  and  periods  and  velocities, 


r (x,y) 


N 

E 

i=l 


(xn:-x)  (yi-y) 


(21) 


where  x and  y are  the  sample  means,  x.  and  yn.  are 
the  sample  standard  deviations. 
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IV.  DISCUSSION 


Miller  (1976)  concluded  that  breaker  shapes  form  a con- 
tinuous gradation  of  shapes  from  spilling  to  plunging.  Each 
breaker  type  undergoes  the  same  processes,  but  differs  only 
in  the  magnitude  of  the  vortices  generated.  The  probability 
density  functions  appear  to  show  a similar  gradation.  Two 
simplified  cases,  for  spilling  and  plunging,  were  chosen  to 
represent  extreme  examples  and  are  used  for  all  discussions. 
The  4 March  1975  data  most  simply  represent  the  plunging 
case  and  the  16  March  1977  data  represent  the  spilling  case. 
The  results  for  all  data  are  given  in  Appendix  C. 

In  both  examples  for  spilling  and  plunging  breakers,  the 
offshore  waves  appeared  as  narrow  band  swell.  A spectral 
analysis  of  the  surface  elevation  records  in  the  surf  zone 
yielded  spectral  width  parameters  ranging  from  0.94-0.97 
(see  Appendix  A),  suggesting  a broad  band  spectrum  in  the 
surf  zone.  The  broad  spectral  width  can  be  attributed  to 
the  non-linear  wave-wave  interactions  yielding  harmonics  in 
the  breaking  zone. 

In  order  to  determine  the  origin  of  the  large  spectral 
width  parameter,  e was  calculated  in  the  ranges  of  0-0.25 
Hz,  0-0.50  Hz,  0-1.0  Hz  and  0-nvquist  frequency.  The  analy- 
sis revealed  that  large  values  of  the  spectral  width  para- 
meter occurred  in  all  cases.  In  the  case  of  0-0.25  Hz  which 
contains  the  first  three  harmonics,  e was  greater  than  0.7, 
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which  was  the  case  with  the  smallest  t . It  is  concluded 
that  the  change  from  a narrow  deep  water  spectrum  to  a 
broad-band  shoaling  spectrum  is  primarily  due  to  generation 
of  secondary  waves  at  harmonic  frequencies. 

A.  PERIOD  PROBABILITY  DENSITY  FUNCTIONS 

Bretschneider  (1959)  made  the  assumption  that  the  wave 
lengths  were  Rayleigh  distributed;  using  this  assumption 
and  the  linear  wave  theory  relationship 


L = -S_ 
2-  1 


(22) 


the  wave  periods  squared  would  be  Rayleigh  distributed  in 
deep  water.  Bretschneider' s analysis  can  be  extended  to 
the  shallow  water  case  using  the  shallow  water  approxima- 
tion for  wave  length 


L = /gK 


(23) 


Assuming  the  wave  lengths  are  Rayleigh  distributed,  the 
periods  in  shallow  water  would  be  expected  to  be  Rayleigh 
dis  tributed. 

The  period  PDFs  for  plunging  and  spilling  examples  are 
shown  in  Figures  5 and  6 respectively.  The  Rayleigh  PDF  has 
been  superimposed  on  the  period  PDFs  for  comparison.  The 
formulas  used  for  the  Rayleigh  PDF  for  various  parameter 
transformations  are  given  in  Appendix  A.  In  both  the  spill- 
ing and  plunging  wave  examples,  it  is  noted  that  the  period 
PDFs  only  loosely  follow  the  Rayleigh  PDF.  The  Rayleigh  PDF 
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overpredicts  the  wave  periods  greater  than  the  Rayleigh  mode 
value,  and  underpredicts  the  observed  periods  less  than  the 
peak  Rayleigh  value.  The  truncation  of  the  long  period 
waves  is  probably  due  to  the  larger  period  waves  breaking 
farther  offshore  and  then  reforming  to  be  recorded  by  the 
stationary  wave  gauge  as  higher  frequency  waves. 

The  spilling  example  for  16  March  1977,  show,  in  Figure  5, 
has  the  mode  lying  close  to  the  maximum  of  the  superimposed 
Rayleigh  PDF.  The  peak  of  the  breaking  wave  distribution 
is  located  at  the  first  harmonic  of  the  offshore  wave  period, 
which  is  attributed  to  the  secondary  waves.  A peak  corre- 
sponding to  the  peak  offshore  wave  period  is  absent. 

The  plunging  waves,  as  represented  by  the  4 March  1975 
record  in  Figure  5,  are  multimodal.  Since  the  Gaussian,  Gram- 
Charlier,  Rayleigh  and  Weibull  PDFs  are  unimodal,  none  of  the 
theoretical  PDFs  will  adequately  describe  multi-modal  PDFs 
for  this  particular  plunging  surf  condition. 

Two  strong  peaks  are  found  to  occur  in  the  plunging 
breaker  period  PDF.  The  longer  period  is  found  to  correlate 
well  with  peak  period  of  the  wave  spectrum.  The  shorter  wave 
period  appears  at  approximately  the  first  harmonic.  Since 
the  center  period  of  the  bin  is  used  to  represent  the  inter- 
val, the  magnitude  of  the  peaks  and  the  proximity  to  har- 
monics as  represented  by  the  labels  is  a function  of  the 
interval  width  and  the  interval  cutoffs.  A 30  bin  PDF  was 
calculated  for  better  resolution  to  determine  if  the  peaks 
were  harmonics  and  Figure  7 shows  that  they  are. 
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ii&ure  7.  PDFs  of  Period  and  Height,  30  bin,  4 March  19  75, 


The  waves  periods  as  determined  by  a zero-up-cross  of 
the  surface  elevation  record  for  spilling  conditions  yields 
a poor  indication  of  the  peak  of  the  energy  spectrum  which 


corresponds  to  the  mean  offshore  wave  period.  In  the  spill- 
ing breaker  case,  the  mean  period  corresponded  to  the  first 
harmonic  of  peak  in  the  energy  spectrum.  For  the  plunging 
breakers,  the  mean  period  corresponded  to  a valley  in  the 
period  PDF  between  the  primary  and  first  harmonic  periods. 

The  appearance  of  strong  peaks  at  one  half  the  offshore  wave 
period  is  indicative  of  the  importance  of  non-linearities 
in  the  form  of  secondar3/  waves  in  the  surf  zone. 

B.  HEIGHT  PROBABILITY  DENSITY  FUNCTIONS 

Wave  height  PDFs  for  the  plunging  and  spilling  wave  ex- 
amples are  shown  plotted  against  the  theoretical  Rayleigh 
PDF  in  Figures  5 and  6.  The  height  PDF  for  spilling  waves 
is  not  as  easily  described  as  the  period  PDF.  The  height 
PDF  is  loosely  described  by  the  Rayleigh  PDF.  The  major 
height  peak  is  centered  at  2.2.  The  Rayleigh  PDF  overesti- 
mates the  density  of  the  small  wave  observations  and  tends 
to  ’underestimate  the  density  of  the  large  waves.  The  heights 
are  truncated  at  approximately  six  standard  deviations.  Since 
the  wave  gauge  is  stationary  and  therefore  the  depth  of  the 
water  is  essentially  fixed,  this  truncation  occurs  due  to 
the  breaking  of  the  large  waves  seaward  of  the  wave  gauge. 

The  calculated  mean  height  value  on  16  March  1977  lies  at 
2.7,  which  actually  corresponds  to  a dip  in  the  PDF  relative 
to  the  peak  intervals  on  either  side. 
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The  mean  height  observed  in  the  4 March  1975  plunging 
breaker  PDF  was  calculated  to  be  at  2.4.  The  observed  PDF 
is  found  to  be  strongly  bimodai  with  peaks  at  1.8  and  3.4. 

The  strong  bimodality  is  attributed  to  the  very  narrow  band 
incident  wave  system  and  the  generation  of  secondary  waves 
at  the  first  harmonic. 

C.  JOINT  PROBABILITY  DENSITY  FUNCTIONS  (T  vs  H) 

The  joint  PDF  of  periods  and  heights  for  the  plunging 
and  spilling  examples  are  contoured  in  Figures  5 and  6;  the 
distributions  for  all  cases  are  found  in  Appendix  C. 

It  is  noted  that  the  longest  wave  periods  are  associated 
with  the  greatest  wave  heights  suggesting  strong  correlation 
between  periods  and  heights.  The  correlation  coefficient 
was  calculated  for  all  cases  and  varied  between  0.60  and  0.80. 
Bretschneider  (1959)  examined  the  two  special  cases  of  joint 
Rayleigh  distributions  with  correlations  of  0 and  +1.  Per- 
fect correlation  between  two  variables  would  result  in 
correlation  coefficient  of  1.0  and  a straight  line  in  the 
joint  PDF. 

The  joint  Rayleigh  PDF  for  the  case  of  a correlation  co- 
efficient of  zero  is  plotted  in  Figure  3.  Zero  correlation, 
or  statistical  independence,  results  in  the  joint  PDF  being 
spread  out  as  compared  with  the  straight  line  for  perfect 
correlation.  The  amount  of  spread  is  a measure  of  the  reduc- 
tion in  correlation. 

Figure  5 is  an  example  of  p(H,T)  for  plunging  breakers 
and  shows  two  regions  of  high  density.  The  bimodality  of 
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wave  period,  and  the  node  of  the  observed  breaker  heights. 

Greater  density  was  found  at  short  periods  and  heights 
for  all  records  than  predicted  by  the  Rayleigh  PDF.  A high 
frequency  wave  riding  on  a long  period  wave  causes  an  addi- 
tional zero  crossing  as  the  longer  wave  passes  through  zero. 
Most  previous  comparisons  of  heights  or  periods  of  waves  have 
been  done  using  pressure  sensors  which  greatly  filter  the 
high  frequencies  due  to  hydrodynamic  attenuation.  Capaci- 
tance wave  staffs  and  electromagnetic  flowmeters  have  good 
response  times  and  thus  measure  high  frequencies.  Therefore, 
even  though  the  data  were  filtered  at  1 IIz , all  records  showed 
peaks  at  the  shortest  periods  and  heights. 
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D.  VELOCITY  PROBABILITY  DENSITY  FUNCTIONS 


The  velocity  PDFs  were  calculated  as  the  maximum  onshore 
velocities  within  the  period  segments  determined  by  the  zero- 
up-crosses  of  the  surface  elevation  record  as  shown  in  Fig- 
ures 8 and  9.  Using  this  technique  can  result  in  both  posi- 
tive and  negative  maximas  for  onshore  particle  velocities 
(i  . e . , d^v/  dt  < 0 ) . 

The  velocity  PDFs  for  spilling  breakers  and  plunging 
breakers  are  shown  in  Figures  8 and  9.  The  theoretical 
Rayleigh  PDF  is  superimposed  on  the  same  graphs  for  compari- 
son. It  is  apparent  that  neither  the  spilling  or  plunging 
case  is  adequately  described  by  the  Rayleigh  PDF.  The  veloc- 
ity PDF  for  spilling  breakers  appears  to  be  symmetrical,  with 
the  mode  value  falling  at  the  bin  with  center  velocity  of 
1.5.  There  also  appears  to  be  a small  peak  value  at  0.3. 

When  compared  with  the  plunging  case,  it  is  found  that  the 
peaks  occur  at  the  same  normalized  velocities.  The  plunging 
wave  PDF  is  bimodal  and  has  relatively  few  occurrences  of 
negative  onshore  velocities. 

E.  JOINT  PROBABILITY  DENSITY  FUNCTIONS  (V  vs  H) 

The  joint  probability  distributions  of  velocities  and 
heights  are  shown  in  Figures  8 and  9.  The  results  show  high 
correlation  between  velocities  and  heights.  The  high  posi- 
tive correlation  suggests  that  the  larger  maximum  onshore 
velocities  are  correlated  with  the  larger  wave  heights. 

Figure  8 shows  the  plunging  breaker  case.  Two  peak  areas 
are  found  to  occur  in  the  joint  PDF  which  reflects  the 
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Figure  9.  PDFs  of  Velocity  and  Height,  4 March  1975. 


bimodality  of  the  individual  PDFs.  The  joint  PDF  for  the 
plunging  waves  suggests  a highly  correlated  case  (0.68)  in 
which  the  largest  onshore  velocities  occur  due  to  waves  of 
greatest  height. 

The  spilling  case  as  shown  in  Figure  9 also  shows  high 
correlation  (0.70).  For  this  example,  there  is  only  one 
major  peak  which  corresponds  to  the  individual  PDFs  which 
are  essentially  unimodal.  The  peak  area  of  the  joint  PDF 
has  been  shifted  toward  higher  waves  and  greater  velocities 
than  the  plunging  case. 

F.  JOINT  PROBABILITY  DENSITY  FUNCTIONS  (T  vs  V) 

Figures  10  and  11  show  joint  PDFs  of  periods  and  veloc- 
ities in  which  the  periods  were  determined  from  the  zero-up- 
crosses  of  the  wave  gauge  record  and  the  coincident  maximum 
onshore  velocities  are  determined  from  the  flowmeter  record. 
The  results  are  similar  to  the  joint  PDF  of  period  and 
heights.  Spilling  waves,  as  contoured  in  Figure  11,  show  one 
main  contour  group.  The  maximum  spilling  probability  is 
found  at  the  intersection  of  period  0.84  and  velocity  1.5. 

The  plunging  waves  as  found  in  Figure  10  reveal  two  contour 
peaks.  One  large  contour  group  is  associated  with  the  point 
representing  the  offshore  period  1.56  and  a velocity  of  1.5. 
The  second  large  contour  is  a function  of  the  first  harmonic 
wave  period  of  0.84  and  a velocity  of  0.3. 

The  negative  values  of  the  maximum  velocity  indicate 


Figure  10.  PDFs  of  Period  and  Velocity,  16  March  1977. 


Figure  11.  PDFs  of  Period  and  Velocity,  4 March  1975. 


The  joinc  PDF  of  wave  height  and  velocity  shows  that  the 
negative  velocities  are  associated  with  the  low  wave  heights 
The  joint  PDF  of  velocity  and  period  shows  that  small  nega- 
tive velocities  are  correlated  with  the  short  periods. 
Therefore,  it  is  concluded  the  negative  velocity  maxima 
are  associated  with  short  periods  and  low  wave  heights.  The 
secondary  waves  occur  after  the  passage  of  the  primary  wave 
crest  and  have  a peak  near  the  trough  of  the  primary  wave; 
the  velocities  associated  with  the  primary  wave  would  be 
negative  maximum  in  this  region  so  that  the  measured  veloc- 
ity, which  is  a sum  of  the  primary  and  secondary  waves,  can 
be  small  or  negative. 
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V.  CONCLUSIONS 


Waves  in  Che  surf  zone  are  a highly  non-linear  process 
which  is  evidenced  by  the  appearance  of  secondary  waves. 

The  secondary  waves  appeared  in  the  PDFs  at  the  first  har- 
monic of  the  primary  period. 

The  joint  probability  density  functions  of  the  breaking 
wave  properties  show  high  correlations.  The  correlation 
factor  ranged  between  0.60-0.80  for  periods  and  heights 
and  suggests  that  the  greater  the  wave  period,  the  higher 
the  breaker  height  will  be.  High  correlation  between  the 
periods  and  velocities  and  the  heights  and  velocities  sug- 
gests that  the  maximum  onshore  velocities  are  a function 
of  wave  periods  and  breaker  heights. 

Breaking  waves  may  be  viewed  as  a gradation  between 
two  extremes,  spilling  and  plunging  breakers.  The  differ- 
ence between  the  spilling  and  plunging  breakers  becomes 
apparent  in  the  joint  PDFs.  The  spilling  breakers  yield  a 
joint  PDF  which  has  only  one  high  density  contour  region, 
whereas  the  plunging  breakers  show  two  high  density  contour 
groups.  The  high  density  contour  groups  appear  to  be  most 
associated  with  the  period  PDF. 

The  joint  PDF  (p(T,H))  of  spilling  breakers  has  a maxi- 
mum density  at  a wave  period  equal  to  one  half  the  peak 
period  of  the  wave  energy  density  spectrum,  i.e.  the  first 
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harmonic.  The  corresponding  velocities  (p(v,T))  or  heights 
(p(H,T))  were  likewise  found  to  be  unimodal  with  a peak  den- 
sity which  is  driven  by  the  first  harmonic  wave  period. 


Plunging  breakers,  at  the  opposite  end  of  the  gradation, 

yield  a binodal  joint  PDF.  The  two  major  contour  groups 

appear  to  be  a result  of  the  bimodal  density  distribution 

of  the  periods.  The  two  peak  densities  of  the  period  PDF 

are  representative  of  the  peak  of  the  wave  density  spectrum 

and  a first  harmonic.  Due  to  the  high  correlation  of  periods 

and  heights , the  corresponding  PDF  of  wave  heights  appears 

to  be  bimodal.  The  joint  PDF  of  period  and  velocity  and  the 

joint  PDF  of  velocity  and  height  suggests  high  correlation 

between  these  variables.  The  maximum  onshore  velocities 
/ 

appear  to  be  correlated  with  both  wave  periods  and  wave 
heights . 

A truncated  P^ayleigh  PDF  was  found  to  loosely  approxi- 
mate the  empirical  period  and  the  height  PDFs  of  spilling 
breakers.  This  truncation  seems  to  be  the  result  of  long 
period  waves  with  their  associated  larger  breaker  heights 
breaking  prior  to  the  arrival  at  the  wave  gauge.  The  mode 
of  the  spilling  wave  period  distributions  corresponds  to  the 
first  harmonic  of  the  peak  of  the  energy  density  spectrum. 

The  smoothed  appearance  of  the  spilling  wave  periods  PDF, 
relative  to  the  plunging  case,  appears  due  to  the  smearing 
which  occurs  with  the  appearance  of  harmonics  greater  than 
the  first.  The  velocities  PDFs  of  spilling  waves,  where  the 
velocity  maximum  is  determined  as  in  Section  D,  yields  a 
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symmetrical  distribution  with  some  negative  onshore  veloc- 
ities present.  Since  the  Rayleigh  distribution  is  undefined 
for  negative  values,  the  resulting  fit  of  observed  data  with 
the  theoretical  Rayleigh  distribution  is  poor. 

The  PDFs  of  the  plunging  breakers  are  poorly  described 
by  either  the  Rayleigh  PDF  or  other  existing  theoretical 
PDFs  such  as  the  Weibull  PDF.  Theoretical  PDFs  are  unimodal, 
whereas  the  empirical  PDFs  under  plunging  conditions  are 
strongly  bimodal.  The  high  correlation  found  in  the  joint 
PDFs  for  wave  heights  and  maximum  velocities  also  results 
in  a bimodal  distribution  both  in  the  velocities  and  the 
heights  PDFs. 

The  mean  wave  period  as  obtained  by  the  zero-up-cross 
method  in  the  surf  zone  is  a poor  indicator  of  the  mean 
offshore  wave  period.  Spectral  analysis  of  the  breaking 
waves  shows  an  energy  peak  corresponding  to  the  primary  off- 
shore wave  period,  but  the  period  PDFs  of  breaking  waves 
have  a strong  peak  at  the  first  harmonic  of  the  wave  spectra 
peak.  Plunging  waves  had  high  density  peaks  in  the  period 
PDF  corresponding  both  to  the  peak  of  the  wave  spectrum  and 
to  its  first  harmonic.  The  result  is  that  the  mean  of  the 
breaking  waves  falls  between  the  primary  and  the  first  har- 
monic wave  period.  The  mean  period  of  breaking  waves  there- 
fore is  a poor  estimation  of  the  mean  period  of  the  incident 
waves . 
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APPENDIX  A 


STATISTICAL  PARAMETERS  CALCULATED 

1,  Due  to  the  variation  of  sampling  frequencies  between  the 
different  experiments,  it  follows  that  some  of  the  other 
statistical  parameters  will  also  vary.  A summary  of  these 
parameters  follows  where  f^  is  the  Nyquist  frequency,  Af 
is  the  resolution.  At  is  the  interval  between  samples  and 
T is  the  length  of  sampled  record  in  minutes. 


DATE 

JLn(Hz) 

Af  (Hz) 

At  (s 

; ample/sec) 

T (min) 

4 

MARCH 

19  75 

2.5 

6. lxlO'4 

0.200 

27.  3 

5 

MARCH 

1975 

2.5 

6. lxlO"4 

0.200 

27.3 

6 

MARCH 

19  75 

2.5 

6. lxlO"4 

0.200 

27.3 

8 

MARCH 

1975 

2.0 

4. 8xl0-4 

0.256 

35.0 

29 

MAY 

19  75 

2.0 

4. 9x10 "4 

0.256 

35.0 

16 

MARCH 

1977 

2.0 

4.9xl0~4 

0.250 

34. 1 

17 

MARCH 

1977 

2.0 

4.9xl0"4 

0.250 

34.1 

19 

MARCH 

1977 

2.0 

4.9x10  u 

0.250 

34.  1 

3 

MARCH 

1973 

2.0 

4. 9xl0~4 

0.250 

34.1 

9 

MARCH 

1978 

2.0 

4.9xl0"4 

0.250 

34.  1 

2. 

Tables  of 

calculated  statistical 

parameters 

follow 

n = 

sea  ; 

surface  elevation 

(m) 

V = 

ins  tantaneous 

particle 

velocities  (m/sec) 

VM  = maximum  onshore  particle  velocities  (m/sec) 

LR  = length  of  record  for  calculated  statistics  (min) 
e = spectral  width  parameter 
H = wave  height  (m) 

T = wave  period  (sec) 
r = correlation  coefficient 
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102.3  102.3  102. 


APPENDIX  B 


MATHEMATICAL  DERIVATION  OF  NORMALIZED  RAYLEIGH  DISTRIBUTIONS 

Wave  heights,  periods  and  velocities  have  been  compared 
with  the  Rayleigh  distribution 

.2 

p(A)  = — * exp  (-  — j) 

4 a 2 


where  A is  the  non-normalized  value.  Since  all  parameters 
are  normalized  in  this  paper,  the  following  conversion  was 
made : 

let  A’  = A/a  , 


then 


,dA 


p(A')  - p (A)  (^i ) A-A'a 
p(A')  = p (A)  = P (A)  a 

for  wave  heights 

p(H)  = — , e:cp(-%(^|-)2) 

4a  n 


,H  s H Q w H ,2, 

p(~)=  — exP(-^C— ) ) 
n n n 


for  maximum  velocities 


P(V)  = — 2 exp(-%(^-)“') 
4a„  v 


,V  \ — v 
P CTv  4cv 


exp  C-%  ) 2 ) 


v 


for  periods 


P(T)  = 


4a" 


exp  (-h(^)2) 
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where 


a = T/2 


exp  (-%( 2^)2) 

exp  (-%  (-)2) 
T 


APPENDIX  C 

PROBABILITY  DENSITY  FUNCTIONS 

The  following  cables  consist  of  both  the  normalized 
joint  probability  density  functions,  and  the  individual 
probability  density  functions  as  labeled.  The  PDF  tables 
coincide  with  Figures  5 thru  11,  respectively. 

The  top  table  is  the  joint  PDF.  The  top  line  and  the 
left  column  represent  the  center  period,  height,  or 
velocity  point  of  each  of  the  15  bins.  The  numbers  in 
the  other  columns  and  lines  represent  the  frequency  of 
occurrence  in  the  respective  bins. 

The  bottom  table  represents  the  individual  PDFs  as 
labeled.  The  top  line  of  the  individual  PDFs  is  Che 
center  point  of  the  individual  bins.  The  second  line 
represents  the  frequency  of  occurrence  in  the  bin. 
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PDFs  of  Period  and  Height,  4 March  1975. 


PDFs  of  Period  and  Height,  16  March  1977. 
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